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DES-TIME-2006 


Dresden International Symposium on 
Technology and its Integration into 
Mathematics Education 

This symposium combines 2 conferences: 

9th ACDCA Summer Academy 

and 

7th Int’l Derive and TI-CAS Conference 

20.-23. July 2006 

Dresden, Saxony, Germany 

1st Announcement and 
Call for Submissions 



Deadline: 15.02.2006 

9th ACDCA Summer Academy 

Topics include: 

- CAS-based curricula and teaching methods 

- Dynamic geometry systems as teaching tools 

- Assessing with technology 

- Internet as a teaching tool 

- Dispensable and indispensable mathematical skills and abilities 

- CAS as pedagogical tools for visualization, experimentation, concentration ... 

- New classroom examples using CAS 

- the next generation: PeCAS=Pedagogical CAS 


7th Derive & TI-CAS Conference 

Topics include: 

- Applications in mathematics, the natural sciences, research & development, eco 
nomy, social sciences, industry, ... 

- Interfaces to other programs / tools 

- Creating and using utility files 

- Programming in DERIVE, in TI-89 / Voyage 200, and in other CAS 

- Problems and limitations 

- Web-Resources - going online 

- Teaching and Learning with DERIVE and other CAS 

- Connecting computer algebra and computational logic 


Please indicate 

- what type of presentation you want to make (lecture, workshop, ...) 

- for which of the two conference strands you intend the submission and submit using the 
respective form at the conference website (www.des-time-2006.de) 
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Dear DUG-Members, 

I am sorry for the delay in delivering DNL#59. It 
was and still is a very busy time for a retired per- 
son. I had the intention to produce a smaller DNL 
this time and make a large issue by the end of this 
year - it will be DNL#60. But as you can see, it are 
again 52 pages. We have two very extended contri- 
butions - one of them from the fields of chemistry. 
I am very grateful to Dr Colin Kennedy for his 
article and especially for his interesting theoretic 
explanation. It's great to have a real life CAS- 
application in this DNL. I could not resist to trans- 
fer ’’Weak Acid Speciation to the TT’. By reproduc- 
ing the functions on the handheld I came across the 
’’Differential Matrix”. Ever heard about this 
(DERIVE function DIF_DATA)? Find a comment on 
it on page 2 1 . 

The second ’’real life application” is modeling a 
Pringle, which seems to be a canned computer 
designed potato chip. Benno Grabinger demon- 
strates how motivated students can be if they are 
facing mathematics in their environment. 

The full German version of Benno 's article can be 
downloaded in the next future from 
www.austromath.at/t3. It is not too difficult to use 
the TIs for this project. 

I received many requests for the User Forum. You 
can find two of them together with the answers 
which might be of interest for you. The remaining 
questions and answers will be published in the next 
DNL. As a reaction on the demand to produce 
contour plots in DERIVE and incited by my at- 
tempts Peter Schofield provided a wonderful pack- 
age which will be presented in the ’’Christmas- 
issue”. 

Furthermore I am proud to announce a contribution 
from Johann Wiesenbauer who taught DERIVE to 
solve SUDOKUs. I must admit that I didn't know 
what SUDOKUs are when he informed me about 
this DERIVE challenge. He is busy providing an 
extremely user friendly program. 

Hubert Weller sent an article about ’’Mathematics 
and Design" and last but not least I received a paper 
from Bernhard Kutzler ’’Numerics versus Sym- 
bolics”. 


I collected a sequence of emails which show how 
CAS-people are thinking and discussing about 
possible developments of a program. Albert Rich 
wrote: 

It brings to light how my thinking evolves in 
response to the dictates of mathematics and 
the comments of real world users like David's 
Stoutemyer and Jeffrey. 

I am very grateful for the permission to publish this 
exciting exchange of mails (pages 40 to 52). Your 
are friendly invited to join the party if you have any 
comments or proposals on this objective. The dis- 
cussion has not come to an end! 

Philip Yorke from England asks for your assistance 
to set up a DERIVE competition (last page). If you 
have any ideas then please contact him directly. It 
would be great to inform me, too. 

Finally I'd like to remind you on our next DERIVE 
Conference DES-TIME 2006 in Dresden. We are 
very proud that we can announce four excellent 
keynote speakers: 

Prof. Kristina Reiss, Augsburg, Germany 

Prof. Willibald Dorfler, Klagenfurt, Austria 

Prof. Carl L. Leinbach, Gettysburg, USA 

Dr. Terence Etchells, Liverpool 

On the first page you can find the First Call for 
Submissions. It would be wonderful to meet you in 
Saxony next summer. Please encourage your col- 
leagues to attend the Conference and to submit 
papers. All contributions are very welcome. 

It remains for me to wish you a fine autumn. I hope 
that you like this issue and I would be happy to 
meet you again in December with DNL#60. 



Download all DNL-Derive- and Tl-files from 

http : //www . austromath . at/ dug/ 
http : //www . bk- teachware . com/ main . asp?session=37 505 9 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS-77 User Group. 
It is published at least four times a year 
with a contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the 77-89 /92/ Titanium! Voyage 200 as 
well as to create a group to discuss the 
possibilities of new methodical and didac- 
tical manners in teaching mathematics. 

As many of the DERIVE Users are also 
using the CAS-77v the DNL tries to com- 
bine the applications of these modem tech- 
nologies. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: December 2005 

Deadline 15 November 2005 


Preview: Contributions waiting to be published (selection) 

Two Stage Least Squares, M. R. Phillips, USA 

Some simulations of Random Experiments, J. Bohm, AUT & L. Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Another Task for End Examination, J. Lechner, AUT 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

ANOVA with DERIVE & Tl, M. R. Phillips, USA 

Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE , J. Bohm, AUT 
Henon & Co, J. Bohm 

Challenges from Fermat, Bj. Felsager, DEN 
Are all Bodies falling equally fast, J. Lechner, AUT 
Modelling Traffic Density, Th. Himmelbauer, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, P. Luke-Rosendahl, GER 

Diophantine Polynomials, Duncan McDougall, CAN 

Rosettes, J. Lechner, AUT 

Jacobi-lteration, Heinz Rainer Geyer, GER 

Numerics versus Symbolics, Bernhard Kutzler, AUT 

Mathematics & Design, Hubert Weller, GER 

Implicit Plots 2D & 3D, Peter Schofield, UK 

and Setif, FRA; Vermeylen, BEL; Leinbach, USA; Koller, AUT; Baumann, GER; 
Keunecke, GER and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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Michel Beaudin, Montreal 

I am writing this email from Nara, Japan (see the Education session of ACA-05 at 

http://www.jssac.org/Conference/ACA/). This is my first trip to Japan: beautiful country, rich history and cul- 
ture! My wife loves it! 

In my talk, I used Derive 6 slider bars for 2D and 3D animation (I should repeat this next year in Dresden). But 
when I used the ’’show step simplification", the audience could not believe it! I used the (known) example of 
int(l/(l+tan(x) A r),x,0,pi/2) whose value is pi/4 for any value of the parameter r. 

I also used the TI Voyage 200 in order to show that, if we want computer algebra to be used by students, we 
must do it each day, in the classroom, not only in computer labs. 

I think that it is very important to show Derive and the TI outside Derive and TI conferences. AC A conferences 
should not be only Maple and Mathematica oriented ... 

Best regards and see you next year in Dresden, Germany! 

Michel Beaudin 


Don Phillips, Washington 

Josef, 

Hello! Hope you've had a good summer!. I've been busy. I put my actuarial math programs from 
Derive into the TI-89. Everything seems to be working as it should. Once I get a write up for how to 
use it, I'll send it on to you. I'm also almost done with putting my bond programs onto the TI-89, cal- 
endar functions and all. I'll send that onto you also when I write up some instructions. 

Regards, 

Don 

...it is unworthy of excellent men to lose hours like slaves 
in the labor of computation. 

Gottfried Wilhelm Leibniz 

Manuel Garcia Mateos, Germany 


2. Preparing a teaching project to the subject „Complex Numbers" I used DNL#55 to among others 
produce a geometric representation of the fundamental theorem. In this DNL I saw that a contour 
plot was made by means of the TI. 

My question is how to project the contour plot into the xy-plane in order to show the graph 
mod_surf(f) in a 3D-Plot Window? It would be great to implement a contour plot command (func- 
tion) as a package in DERIVE (similar to MAPLE). 

Many thanks 
Manuel Garcia Mateos. 


DN.: Manuel sent a great article - translated into English - which I forwarded to Bert Waits and Frank 
Demana, who were the authors of the contribution in DNL#55. This article will be published in one of 
the next DNLs , many thanks Garcia. Bert immediately wrote back: 


Many, many thanks for 1 . Publishing the paper and 2. sending this email! It is great, as you noted, 
to see ones work recognized so nicely! 

Bert 


p 
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But this was not the whole reaction. I sent the request for contour plots to Albert Rich and Theresa 
Shelby who involved our friend Peter Schofield. 

In the meanwhile the problem how to provide contour plots made me busy and after some attempts I 
found a way to produce them and sent my first results to Garcia, Albert and Peter (see the plots on the 
bottom of this page). This is Albert's answer: 

Very clever. Looks like it has the makings of a nice article for the Newsletter. 

Peter studied my function and it made him busy, too: 


Hello Josef, Albert and Theresa, 

Thank you for your email Josef. 

Your contour plotting idea is extremely interesting. It reminds me of the glass ornament cubes (you 
see in shops) with the lazar etched figures inside. The finished 3D-plots are very attractive. 

The calculations appear to take a long time. I am attaching two Derive files that may suggest ways of 
speeding them up. 

1. polycontour.mth 

Calculation times can sometimes be reduced by restricting your function input to implicit polynomial 
equations of form f(x,y,z) = 0. This file also contains another way of persuading Derive to 3D-plot 
small points. (It still takes the 3D-plot widow around 60 seconds to process some of the large data 
lists.) 

2. Nsolutions.mth 


This file contains a faster procedure (N_SOLUTIONS) for finding all of the zeros of f(x) = 0 
for a < x < b. I’m not entirely happy with it, but it works. Perhaps Albert can improve it? 


In fact, something like this for both NSOLUTIONS and NSOLVE would be a useful addition to a futu- 
re update of Derive 6. 


Can I request it for the wish list? 
Hope this helps. 

Best wishes, 


Here are two results of contour plots. Another one is 
given at the end of the “Pringles” contribution. 

We will come back with the full contribution in the 
next DNL! 


Peter 
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Dirk Callens, Belgium 

Hi Josef, 

I have a little problem with Derive. (See attachment). Can you solve it ? 
Thanks, Dirk 


#1: CaseMode := Sensitive 

We have two points A and B, and a vector V of points in the plane. 

1 l 

2 2 

3 3 

4 4 

3 3 


# 2 : 


A := [4, 2], B := [8, 2], V := 



a) Sort the vector V in descending ordre of the distance to the origine. 

3 5 


#3: SORT(V, |X| > |Y| , X, Y) = 


Ok.this works. 


4 4 
3 3 
2 2 
1 1 


b) Sort the vector V in descending ordre of the distance to the point [2,2] and 

c) Sort the vector V in descending ordre of the surface of the triangle with vertices A, B and a point of V. 


#4: SORTCV, |X - [2, 2]| > |Y - [2, 2] | , X, Y) = SORT 


r 

"11" 



2 2 



3 3 

, ? , X, Y 


4 4 



.33. 

J 


#3: SORT 


1 1 

V, |CRQSS(X - A, X - B) | > |CRQSS(Y - A, Y - B) | , X, Y 

2 2 


SORT 



"11" 



2 2 



UJ 

LU 

, ? , X, Y 


4 4 



.33. 

J 


Why don't both procedures work ?? 

Although we can make a table of the areas ! ! 


#6: TABLE 


-■ |CR0SS(X - A, X - B) | , X, V 


[ 1 , 1 ] 2 
[ 2 , 2 ] 0 
[3, 3] 2 

[4, 4] 4 

[5. 5] & 
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I sent Dirk's request to Johann Wiesenbauer, our “Titbits Man” and he very quick found a satisfying 
solution for sorting under conditions: 


Conditional Sorting with Derive 


As for conditional sorting, it is rather the rule than the exception that Derive can’t cope 
with user defined comparisons. 

It can though, with "a little help” of the user. The trick is to insert the expression used 

for the sorting as a first component and to sort this two element vector. Then the second column 

gives the result. This is what my function mysort(v,a,e) does, where v denotes the vector 

to be sorted, e is the expression that is used for sorting and a stands for sorting "in ascending order". 


mysortCv, e := x, a := true) := 

If a 

#1: (SORT (VECTOR ( [SUBST (e , x, v_), v_], v_, v))) COL 2 

REVERSE (myso rt(v , e) ) 


f S 5 1 


#2: mysort(v, |x - [2, 2]|, false) 


4 4 
3 3 

1 1 


L 2 2 J 

r 2 2 1 


#3: 


mysort 


1 

|CR0SS(x - A, x - B) 

2 


1 1 

3 3 

4 4 


L S S J 


f 5 5 1 


#4: 


mysort 


v, 


1 


|CR0SS(x - A, x - B) 

2 


false 


4 4 
3 3 

1 1 


L 2 2 J 


Shaded Riemann Sums 

Martin Wuttke from Berlin wrote to Walter Wegscheider, the author of the DERIVE Online Tutorial 
(www.austromath.at/daten/derive) how to produce the shaded area under the step-function (Riemann 
sums). He tried working with AreaUnderCurve and found it very tiresome doing it piece by piece. Wal- 
ter had produced the illustration in his tutorial using GeoGebra (www.geogebra.at) and asked me if I 
had any easier way to do this with DERIVE 

DNL: See the results of my efforts on the next page. You will find how to plot "real" upper- and lower 
sums in one of the next revised old DNLs. 
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2 

#22: osumsh(x , -3, 4, 14) 



#5: rsuni5h(SIN(x) , 0, tt, 10) 


#6: SIN(x) 



2 

#8: msumshCx , 0, 4, 4) 



The Functions above are small programs, the next function is one VECTOR-command in- 
volving the CHI-function: 


msumchi' (u , a, b, n) : = 


VECTOR 0 < y < SUBST u , x , a + 


b - 


(b - a)- (k - 1) 


(b - a) - (k - 1) 


k- (b - a) 


X * + 


k, 1, 


u 


msumchi (SIN(x) , 0, n, 20) 



ri emsum. dfw 
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Download the new Operating System for TI-CAS devices! 


Voyage200/TI-89T OS enhancements 

Operating System Version 3.10 


Key New Functionality 

• Graphing Improvements ► 

• Domain Improvements 

• Gradian Angle Measure 

• Population Standard Deviation 

• Implicit derivatives 

• Increased Solving Capabilities 

• Log to Any Base and Root Functions 

• Bugfixes 

Availability 

• Available next from the Tl web site 

• Available “out of the box”: Running change, not available pre- 
loaded on units approx, until January 1 st TI-89T, April 1 st Voyage 
200 

Compatibility 

• Compatible with Voyage 200 and TI-89 Titanium! 


Too Is ZOOM 




^yl=tan( x) 
y2= 

„ *h. 



I ! 








■[-6 0]4a [-6 0] 

■[-1 2] -»b [-12] 

■ jolw(x-a + y-b=[7 6],C*> 
x = - 5/"3 and y = 3 

■s-olv»g(x 2 < l,x) 




logjClOD- log 3 (5) 


lo-g 3 <2) 


Customer needs addressed 

- New OS extends current CAS functionality to better meet advocates’ 
expectations, especially those which our competitors have already 
addressed. 

- New OS corrects some of the long-standing deficiencies of the CAS. 

Key Benefits 

- More Pedagogical: Domain Tracking, Graphing and other 
Improvements make our solution a better then before fit for high 
school math environment. 

- More Powerful: Solving Inequalities, Implicit Derivatives and Gradian 
Angle Measure make our CAS system more versatile. 

- Easier to Use: Ease of use features include easier access to Apps 
from Desktop, better interface for solving systems of equations using 
vectors, root notation and log to any base functionality. 

- Support for European Audience: The new OS is supported by 
updated localizer Apps and comprehensive guidebooks. 


• Block auto-simplification when 
graphing 

• Warn the user that intermediate 
steps in a computation involved 
complex numbers and when the 
domain may have changed 

• Give the user the option to 
remove “faux” asymptotes by 
introducing a smart discontinuity 
detection mechanism. 

• Graphing piecewise defined 
functions improved as well! 
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• Support for logarithm to any base 
including transfinite bases and 
algebraic expressions. 

• Easier access to log functionality. 

• Try to prevent the user noticeable 
conversion to natural log. 

• Natural log and any log conversion 
commands. 

• Subscripts. 

• Support for root function 


r.‘fflRl 9 r Sr.IcSrcl0t r KrlPr5,I0lcu3T UpI 1 

■ log x (0) 

undef 

■ log 3 C10) - log 3 (5) 

log 3 (2) 

■ log 5 (x + 2) + log 5 (x) | x > 0 

log 5 (x (x + 2)) 

■ ( log 3 (2))> logbase( 10) 

log(2) 

log(3) 

losr<2.3> ►lo«rbase<10>l 

MAIN RAO AUTO FUNC S/30 


| ▼ |r 1 gebr alca^l c]other]pr gn I ojc 1 ean 

Up| 1 


1 

X* 

- 3 J3 1, 

.44225 

■ 3 J8 

2 

■ 3 J~~8 

-2 


MAIN GRO AUTO FUNC S/30 


• Increase the ability of the Tl CAS 
to cover inequalities involving 
polynomial, rational, and absolute 
value expressions. 

• Solve fundamental problems 
involving vectors. These include 
decomposition and intersection 
problems (two lines, a line and a 
plane, two planes) involving 
parametric forms. 


(w ^iRlJIbrJc^clot^rfpr^Iolc 1 ean UplPl 

■ soWeCx 2 > L , x) M or x>l 

■ soluej^- > 1 f xj x#Qand’l<x<l 

■ solt»#(x 2 > * 1 , x) true 

main mi. auto fuhc a/to 


b* f— |r l gebr a |ca 1 clo therjpr gn I ojc 1 ean Up! 1 

■[-6 0 1] ■> d 

[-6 0 1] 

■[4 O 2] •> a 

[4 O 2] 

■ [-1 2 l]*b 

C-l 2 1] 

•17 6 3] + c 

17 6 5] 

■ solve(x a ♦ y b + z-c 1 

= d, Cx y z» 

x = 11/6 and y = 4 and z = * 4/3 


MAIN RAO AUTO 

FUNC 5/30 


• Gradian angle measure has 
long been a request from 
European teachers. 

• Add Gradian angle MODE. 

• Add Gradian angle indicator on 
the status line: GRD. 

• Use superscript G to indicate 
the new angle measure. 

• Implement ►GRAD for angle 
conversion. 



MAIN RAO AUTO FliNC 3/30 


i^Ifi 1 gebr alca rdotherlprgn I ojc 1 ean Up! ] 

■ sir* 400 G ) 

0 

■(345° )>6rad 

(-L^G 

■(345° )A6rad 

383. 333 Q 

■ ((2/3 • n) r )>Grad 

m 

■((2/3 Ji) r > Grad 

133. 333 g 

< <2/3*n> r > ►GradJ 

MAIN GRD AUTO 

FUNC 5/30 
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d{C SB ) _ 

d (K) (v f +v b ) 2 


Writing Vf as V h the volume added values in the first column of the matrices FDMP, FDDP and 
FDTP, then the buffer capacity can be written as- 

„ d{C, A die,,) dv : C R V f dV. 


djCss) _ djCsn) dV t _ C B V f dV^ 
d{pH i ) dV t d(pH t ) (v f+Vi ) 2 d(pH i ) 


So, p = 



1 

f d(pH) 

) 


The Derive function DIF_DATA has already been used to calculate a 2 column matrix in which the 

( d (pH 

first column is V,- values and the second column, 1 1 « k ~ 


. Hence the buffer capacity graph can be 


obtained via (28). The P values are orders of magnitude too small to be plotted on the pH titration 
graph so they are first scaled by a factor ‘multbuf . The buffer capacity curves can be created and 
scaled using the utility file functions BUFCAP_MP(lencol, multbuf), (#52), 
BUFCAP_DP(lencol, multbuf) (#56) and BUFCAP_TP(lencol, multbuf) (#60) for mono-, di- and 
triprotic acid titrations respectively. 

Reference 

1. Pankow, J. F., 1991, Aquatic Chemistry Concepts, Lewis Publishers, Chelsea, Michigan, USA, pp 46-49. 


(2) Weak Acid Utility File: Creating Speciation and Titration Graphs for 
Mono-, Di- and Tripodic Acids 

File Weak Acid utility file CDK.mth is among the files connected with this DNL and can be 
downloaded. 


(3) Some Examples using the Utility File functions: 

Malonic Acid, a weak diprotic acid. 

#1: L0AD(C: \ Prog ram Files\TI Education\Derive 6\Users\\Weak Acid utility file 

CDK.mth) 

It is supposed that we have a 0.002 M solution of malonic acid which also contains the 
strong electrolyte KCI at a molar concentration of 0.1 . Except at extreme low and high pH, 
the ionic strength, I, will be determined by the KCI present. So- 

#2: cone := 0.1 






D 


-N-L#59 


Colin Kennedy: Weak Acid Speciation a. o. with Derive 


p!7 


#3: num_cat := 1 

#4: cat_charge := 1 

#5: an_charge := -1 

#6: SING_SALT_I(conc, num_cat, cat_charge, an_charge) 

#7: 0.1 

#8: I := 0.1 

The 'thermodynamic' dissociation constants for malonic acid at 298 K are- 
#9: pKtl := 2.85 

#10: pKt2 := 5.7 

At 298 K, A = 0.51 , the charge on the fully protonated acid, 'z' is 0, and for the first dissociation step, n 
=1 . So the first 'mixed' acid dissociation constant is given by- 

#11: A := 0.51 

#12: PKMNCpKtl, A, I, 0, 1) 

#13: 2.737670932 

The second 'mixed' acid dissociation constant is- 
#14: PKMN(pKt2 , A, I, 0, 2) 

#15: 5.363012797 

So, to 3 significant figures- 
#16: pKml := 2.74 

#17: pKm2 := 5.36 


Molar Concentration versus pH Graphs 

For a 0.002 M solution of the acid the individual species concentrations may now be plotted as a func- 
tion of pH. (Black curves in the graph) 


#18: 

CA := 0. 

002 


#19: 

H2A(CA, 

pKml, 

pKm2) 

#20: 

HAm(CA, 

pKml, 

pKm2) 

#21: 

Am2(CA, 

pKml, 

pKm2) 


At very low ionic strength, the pKm and pKt values will be almost equal, so the speciation curves would 
then obtained from- 


#22: 

H2A(CA, 

pKtl , 

pKt2) 

#23: 

HAm(CA, 

pKtl , 

pKt2) 

#24: 

Am2(CA, 

pKtl, 

pKt2) 


These are the red curves in the graph. Note that a change in ionic strength can make a significant 
difference to the species concentrations. 
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pH Titration Graphs 

Suppose that (Vf = 20) ml of (CAO = 0.002) M malonic acid are titrated against (CB = 0.02) M 
NaOH solution in the burette. All solutions contain 0.1 M KCI, so the ionic strength I, is still 
0.1 . The activity coefficient of hydrogen ions, yH, at this ionic strength is first required- 


#25: 

ACT_COEF_H(A, I) 


#26: 


0.7720953423 

#27: 

yH := 0.772 



The titration curve can now be calculated as a two column matrix (The black curve on the 


graph) 



#28: 

Vf : 

= 20 

#29: 

CAO 

:= 0.002 

#30: 

CB : 

= 0.02 

#31: 

pKw 

:= 14 

#32: 

VOL. 

_BUR_DP(Vf , CAO, CB, pKw, pKml, pKm2, yH) 


The first differential matrix of this is calculated (FDDP) and the column length, (lencol) 
found. - 

#33: FDDP 

#34: lencol := DIM(FDDP COL 2) 

#35: lencol := 101 
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The vertical height of the first differential plot may need to be scaled to fit it on the same 
graph. Possibly a separate scaling may be necessary for each endpoint (red curves). 


#36: multdif := 0.1 

#37: FDDP_SCALED(1 encol , multdif) 


Changing multdif to 6 gives the higher red curve, revealing the first less clear endpoint. 

The buffer capacity of the solution in the titration flask changes throughout the titration, be- 
coming higher as free hydrogen or hydroxyl ion concentrations increase, and often reaching 
a peak at or near half equivalence point titres (as in commonly used buffer solutions). This 
has to be scaled by a large factor (multbuf) to make it visible on the graph (blue curve). 

#38: multbuf := 5000 

#39: BUFCAP_DP(1 encol , multbuf) 



pH titrations are often used to find the pKm's and hence the Km's of weak acids by reading 
off the pH at half equivalence points, here at 1.00 ml and 3.00 ml of titrant added. This 
method yields apparent values read from the graph of pKml = 3.35 and pKm2 = 5.38, in- 
stead of the actual values of 2.74 and 5.36. So the estimation of the second dissociation con- 
stant exponent is reasonable but the first is significantly in error. Discussion with students as 
to why the method fails helps to encourage (1), better understanding of the solution chemis- 
try involved and (2), a more judicious approach to experimental methods. 
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Log Concentration versus pH Graphs 


It is often useful to show speciation graphs as plots of Iog10(molar concentration) versus pH. 
One advantage of this approach is that the full range of hydrogen ion and hydroxyl ion molar 
concentrations can be shown on the same graph as that of the other species present. The 
functions to be plotted are- 

For [H+] (red line) 

#40: -pH - LOGCyH, 10) 

For [OH-] (blue line) 

#41: -pKw + pH - LOGCyH, 10) 

For the 3 species derived (pun not intended!) from malonic acid (black lines on the graph) 


#42: 

L0G(H2A(CA, 

pKml, 

pKm2) , 

10) 

#43: 

L0G(HAm(CA, 

pKml, 

pKm2) , 

10) 

#44: 

L0G(Am2(CA, 

pKml, 

pKm2) , 

10) 
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The concentration values below pH 2 and above pH 12 are less accurate as the ionic 
strength of the solution is no longer largely determined by the 0.1 M KCI in the solution. Nev- 
ertheless, relationships between [H + ] and [OH ] with the concentrations of the species de- 
rived from malonic acid are made apparent from these graphs. 


As my main intention is to show that Derive can be a useful tool in teaching some aspects of 
solution chemistry, I have provided only a few practical examples. However, there are an 
almost infinite number of ways in which the generation of graphs similar to those above can 
simulate practical experiments as well as stimulate an understanding of the principles in- 
volved. 
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Comments on the “Differential Matrix” 

The author of the “Weak Acid ...’’-contribution mentions a “Differential Matrix”(#33). Inspecting his Utility file 
we find that he applies the DERIVE function DIF DATA which is part of NumericalApproximation.mth. 
This is a nice way for numerical differentiation if a function is given by a consecutive sequence of data points. 
We can read in the Online - Help: 

DIF_DATA(A) approximates to the first derivative of a numerical data matrix A. A must have two col- 
umns and at least three rows, with each successive row being a successive [x, y] coordinate pair. 
The result is an analogous data matrix wherein the second column is the approximate first derivative 
of y with respect to x at successive values of x. The x values can be unequally spaced. They can 
even be nonmonotonic, such as samples from a parametrically-defined curve. The formula is based 
on the derivatives of parabolas that fit three successive points. 

It might be useful for readers (and for students) to doublecheck the procedure: We take three consecutive data 
points and find the parabola passing through the points. The 1 st derivative in the 2 nd point is an approximation 
for the derivative of the unknown function containing all data points. There are two exceptions: For the first 
point one has to take the derivative of the parabola given by the first three points in the first point - and similar 
for the last point. Let's check with DERIVE: 


We create the matrix of data points (101 # 86 : 

points) and name it m. : 


The quadratic regression line through the 
first three points serves for two purposes: 

#89: 

It gives the slope for the first and for the 
second point. 

#90: 

#91: 

Points #2, #3 and #4 define a parabola and 

#92 : 

its first derivative in the middle point is 
the slope of the unknown function in point 
#3. m3: 

DIF_DATA uses another algorithm but # 94 ; 
comparing the first lines of DIF_DATA(m) 
which is FDDP in the contribution shows 

#95: 

the correctness of the procedure. 

#96: 


m := V0L_BUR_DP (20 , 0.002, 0.02, 14, pKml, pKm2, 0.772) 
m 

' -7.674608567 2 ' 

-6.547594726 2.1 

-5 .483814654 2.2 

-4.498347975 2.3 


pi := FIT Lx, a-x - 1 - b-x + c_|, |"m , m , m 1 

L L 1 2 3 JJ 

2 

pi := 0. 002407500835 -x + 0. 1229700212 -x + 2.801945902 
d 

— pi = 0.004815001670.x + 0.1229700211 
dx 

SUBSt[ — pi, x, m I = 0.08601676813 
Ldx l,lj 

SUBSt[ — pi, x, m I = 0.09144334165 
Ldx 2 ,lj 


p3 FIT Lx, a-x + b-x + cj, [m , m , m 1 
L L 2 3 4 JJ 

fd ^ 

SUB ST — p3, x, m = 0.09788232195 
3 , 

DIF_DATA(m) 


See the Tl-realisation on page 39. 


-7 .67460856 7 0. 086 01676815 

-6 .547594726 0 . 09144334167 


-5.483814654 


0.09788232195 
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Wie war's mit einem Pringle? 

A Pringle is a Mathematical Potato Chip! 

Benno Grabinger, Neustadt/Weinstralie, www.bennoqrabinqer.de 



Wie kann die Form eines Pringle mathematisch beschrieben werden? 
Wo entstehen solche Fragen? 

Eine Kursfahrt nach Florenz und Siena war der Ausgangspunkt: 


How can we describe the form of a Pringle mathematically? 
Where are such questions coming from? 

A class excursion to Florence and Siena was the initial point. 
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Gipsmodelle (hergestellt von 
Sebastian Schulz, Leibniz- 
Gymnasium Neustadt, 2004) 
erleichtern das Arbeiten. 

Gypsum models (produced by 
Sebastian Schulz) support 
imagination. 




Offenbar sind zwei Symmetrieebenen vorhanden, welche das Pringle langs der gestrichelten 
Linien schneiden. 

Der Einfachheit halber soil einmal angenommen werden, dass die eingezeichneten Schnitte 
der Symmetrieebenen mit dem Pringle Parabeln sind. Wahlt man den Schnittpunkt der 
eingezeichneten Linien als Ursprung eines geeigneten Koordinatensystems, so konnen die 
Gleichungen der Parabeln ermittelt werden, wenn man das Pringle ausmisst. 

Obviously there are two planes of symmetry, which intersect the Pringle along the dashed 
lines. 

To make things easier we will assume that these dashed intersection lines are parabolas. 
Choosing their intersection point as origin of an appropriate system of coordinates we should 
be able to find out the equations of the parabolas by measuring the Pringle. 

(For the following I'll try to translate Benno' s paper into English to save space. Josef. As I 
wrote in the Letter of the Editor, the full German version can be downloaded.) 
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We take the bearings of the center of symme- 
try of the pringle and we obtain point R' 

(0/3/1 .5) as projection of the center on the 
y-z-plane (related to the x'y'z' system of 
coordinates in the vertex of space). 

Points P and Q' are P(2. 3/0/2) and 
Q'(4. 5/3/0). 



As we want to work with a system of coordinates xyz with its center in the center of symmetry 

f Z3 l 

of the pringle and axes parallel to x'y'z' we have to consider the translation v = 3 . Then 


v 1.5 y 


points P and Q have coordinates P(0/-3/0.5) and Q(2.2/0/-1 .5). 


The parabola passing point P lying in yz-plane has the general form z = ay 2 and from 
0.5 = a * 9 follows that a = 1/18. 


A 

So finally we obtain z = jqY 2 


or in parameter form x = 


0 

t 

.18 , 


Very similar we find the parabola with vertex Q lying in the xz-plane, z = b x 2 : 

-1 .5 = b * 2.2 2 which leads to b = -75/242. 

r \ 


z = _25_ x 2^ z = L_ 

242 3.227 


x 2 — > x = 


s 

0 

s 2 

3.227, 


Plot of the parameter representations. 
Derive file modell.dfw 
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Which is the next step? 


Inspecting the intersection curves of the pringle with planes parallel to the yz-plane one get 
the impression that these curves might by also parabolas and in the easiest case these pa- 
rabolas could be the same as the parabola in the yz-plane. In order to reach this goal in the 
model one has to shift this first parabola in such a way that its vertex glides along the parab- 
ola lying in the xz-plane. 


This is achieved by adding a translation vector to the parameter form of the parabola being 
shifted. We find the translation vector by reasoning how the origin is moved. As the origin 
shall glide along the second parabola we have: 


"\ 


( \ 


( 

A 


( 

\ 

s 


0 



S 



s 

0 

and further x = 

t 

+ 


0 

= 


t 

S2 


f2 



S 2 


f2 

s 2 

V 3.227 J 


CO 


V 

3.227J 


ll8 

3.227J 



You can see a family of parabolas which are translated with an increment of s = 0.2. 
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This surface is very similar to the pringle but the boundary lines are missing: 

Pringles must fit into a can in the shape of a cylinder, so we have the idea to intersect the 


surface z = 

18 242 


with a cylinder which has the z-axis as its axis. 



Observing the intersection of cylinder and surface we obtain the Pringle: 



(Expression #6 shows how to "produce" the Pringle. For this purpose we use the rectangular 
form of the cylinder. We don't show the mesh lines and choose an appropriate color in order 
to plot a realistic "naturally coloured" picture of the Pringle.) 

y 2 75 X 2 y2 y2 

The equation which is the base of the Pringle z = - — — = f— - - —— is a special case 

H a 18 242 18 3.227 K 

y 2 x 2 

of the so called hyperbolic paraboloid z = 

One can study the influence of constants a and b for the form of this surface by introducing 
slider bars. 


Derive file hyperbolic paraboloid.dfw 
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File Edit Insert Set Options Window Help 


M 


D^Qi x A S’ 'i 


T i 



Switch to the 3D-Plot Window and issue Insert > Slider bar. 
Fill in as shown below: 


Slider Bar Properties 


Variable Name: 

l a 


Minimum Value: 

F~ 


Maximum Value: 

F~ 


Intervals: 

F" 


W lUpdate plot while sliding! 


OK 

Cancel 

1 

1 1 


Repeat this step and introduce a second slider bar for b. 

Then plot the following expression and use the slider bars to vary the 
values for a and b. 


# 1 : 



1.00 


! 6.00 

J 

□ 

1.00 J 

6.00 





5 Eye: 16; 24, 12 


Center: 0,0,0 


Length: 4.6 : 6 : 6 


McAfee-Produkt-Setup DNL59.doc - Microsof, 


■ ^ Adobe Photoshop Ele. . . 


The Hyperbolic Paraboloid 

The word "paraboloid" is self explanatory because of the gliding parabolas which form the 
surface. The attached adjective "hyperbolic" is based on the form of the intersection curves 
of the surface with horizontal planes. 



It seems to be very likely that the intersection curve is a hyperbola. Taking the plane z = k 

V 2 x ^ 

then the intersection curve is given by k = 



p28 


Benno Grabinger: A Mathematical Potato Chip? 


D-N-L#59 


Examples 


1 . For k = - 1 (picture above) we find: 

,2 


’ ■’ b 2 ' 


a 2 b 2 b 2 a 2 

And this is the equation of a hyperbola in the xy-system of coordinates. 

2. For k = 1 (left graph in the next picture ) we find: 

1 = - — which is again a hyperbola in the yx-system of coordinates. 

3. For k = 0 (right graph in the next picture) we have a special case: 

°=£±£. 

a b 

i.e. the intersection curve falls into two lines. 



These two lines are called the main vertex generating lines ( Hauptscheitelerzeugende ) of the 
hyperbolic paraboloid. Their intersection point lies on the axis of the surface and is called 

vertex. 


The following properties of this surface can immediately be derived from its equation 

_ y 2 x 2 
Z “ a 2 b 2 ' 

There is exact one point of the surface above each point of the top-view-plane (xy-plane). 

The yz-plane x = 0 is symmetry plane of the surface, because variable x occurs only in quad- 
ratic form, i.e. if point P(xlylz) is point of the surface then point P*(-xlylz) is also a point of the 
surface. 

The xz-plane y = 0 is symmetry plane of the surface, too, because variable y occurs only in 
quadratic form, i.e. if point Q(x/y/z) is point of the surface then point Q*(x/-y/z) is also a point 
of the surface. 

The intersection line of the two planes of symmetry - in our presentation the z-axis - is an 
axis of symmetry of the surface. It is called Axis of the hyperbolic paraboloid. 




Derive file intersection. dfw 


There are a lot of applications of this surface in architecture. 
See some examples: 



Saddle Surface: Membrane-roof for the Bundesgartenschau in Cologne, 1957 
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On the website http://www.stellenboschwriters.com/ one can find: 

What is a Hyperbolic Paraboloid? 

Very easy, it is the mathematical name for the form of a saddle. 
Take a piece of textile. 



Bend vertices A and C upwards and fix them. Then bend vertices B and D downwards. If 
your textile is large enough you can use ist as a sun sail as we did on board of our yacht. 
This sun sail shows the same form. It is called HP-surface (HP-Schale). 



Another example: 

http://www.binqham.nl/en/enfotozadel.html 

Saddle Surface 

The saddle surface consists of two maximum- and two mimimum points with a membrane 
tightened between them. 



Zadelvlak ©Tentech 
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The given examples deliver another approach to the hyperbolic paraboloid. We will show that 
the form produced by the described method of hangig up a textile gives a paraboloid. 


y2 

Take the paraboloid z = fy-yy- Four points A, B, C and D on the surface are selected in 

such a way that A and C have z-coordinate -h and y-coordinate 0. Points B and D are above 
the xy-plane with z-coordinate h and x-coordinate 0. 

(See the graph below with hidden point C 


z 



B is defined by B(0/ylh). 

v 2 

For calculating y we substitute the coordinates of B in z = yy 


x^_. 
k 2 ' 


y 2 = hk 2 or \y\ = k4h. So B(0/kJhlh) and D(0/- kjh I h). 


In a similar way we obtain f\(k4h 101- h) and C(-k4h 101- h). 


If we intersect the paraboloid with a cube defined by its vertices A, B, C and D then the re- 
sulting form is very similar to the suspended membrane which we mentioned earlier. 



How to obtain this intersection figure? 
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Point E has coordinates E(0//cV/7 /- h). 

The plane defined by A, B and E: r =r A + AB + AE . 


y 


(kjh) 

0 

+ 171 

1 

V 

+ n 

f-kS] 

k4h 

w 


-h 

\ J 


2/7 

V J 


0 

V J 


This leads to x + y = k\lh. In the same way one can calculate the equations of the remain- 
ing planes of the cube. We collect the conditions for the inner part of the cube: 


|x + y\ < k\[ti and \x -y\< k4h and |z| < h 
Derive file suspended.dfw 


The last figure leads to the conjecture that lines g(A,B) and g(C,D) are lying on the surface - 
that they are part of the surface. This conjecture is checked by calculation: 


g(AB): 

y+_>+ 

k 2 k 2 



(kjh) 


f-kS] 

r = 

0 

+ 171 

k\[ti 


-h 


2/7 


\ J 


b J 


V 2 

Substituting the coordinates into z = 


mWh _ (k4h^k]hf = m2h _ _ m4Ef = _ h + 2mh , 


X 2 

— results in: 


i.e. all points of 


line g(A,B) are elements of the hyperbolic paraboloid. 

Similar calculation shows that all points of line g{C,D) are also elements of the surface. The 
result of the next consideration is, that there are more straight lines which are elements of 
the surface. 

1 . Set up the equations for segments AB and DC. 

2. Subdivide AB and CD into 10 segments of equal length. 

3. Set up the equation for segment PQ, with P being the i th division point of AB and Q 
being the i th division point of DC. 

Then line g{PQ) is part of the surface. 



A 


D 
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Perform this procedure for all division points. 



The proof: 


Segment AB: x = 

fkS] 

0 

+ s 

f-kS] 

k\[h 


-h 

\ J 


2/7 

^ J 


with s e [0; 1] , divisionpoints s = 0.1 , 0.2, ..., 0.9 



f 0 ^ 


f-kyfh) 

Segment DC: x = 

-kih 

h 

V J 

+ s 

ki 

-2/7 

V J 


with s e [0; 1] , divisionpoints s = 0.1 , 0.2 0.9 

Straight line connecting division points P s and Q s (with s fixed and t variable) 


X = 

(k4h) 

0 

+ s 

f-ks] 

kit 7 

+ t 


f 0 ^ 

-k\[h 

+ s 

f -k) 
kit 7 


[kih\ 

0 

-s 

(-kih Y 
kit 7 


—h 

V J 


2/7 

V J 



h 

\ J 


- 2/7 

V J 


—h 

V J 


2/7 

V )_ 


(kjh) 

0 

+ s 

f-kS] 

ki 

+ t 

f-kjh) 

-kiti 

+ st 

f 0 ^ 

0 

-h 

J 


2/7 

V J 


2/7 

V J 


K -Ah y 


This set of points is completely part of the HP-surface because 

x = ki -skih -tk\[h 
y = skih -tk\[h 
z = -h + 2s h + 2th -As th 

y 2 x 2 _ k 2 (s\[h -t\fh ) 2 k 2 (i) -S\fh -t\lh ) 2 _ 

~k 2 ~~k 2 ~ k 2 k 2 " 

= h(s 2 -2st + t 2 -(( / \-s)-t) 2 ) = -h-2sh + 2th + st (-Ah) = z 


What happens if keeping t fixed and s variable with s e [0;1 ]? 
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Answer: for one fixed t 


for various values for f: 




Varying s and t in the interval [0;1] results in 
two families of lines lying completely in the 
surface. 


Better than a graph is a model, of course. 
(This model was built by Dirk Barthel, a stu- 
dent of the Leibniz-Gymnasium Neustadt, 
2004) 
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Changing the discrete steps for s and t to a continuous transition we find the graph as follows 
in the left graph below. 

w 2 y2 

It cannot be distinguished between this surface and z = — - (* * *) which can be 

shown by adding its graph to the existing one (right graph below). 



Derive file membranl.dfw 



The plots are illustrating that both sets of points defined by (**) and (***) are identical. 

The proof by calculation is very easy: 

We showed above that each point of (**) is also point of (***). We will add that each point of 
(***) is also point of (**). 


Let P 


/ / / 

xlyl ^~w 


2 \ 


an arbitrary point of (***). Parameter form of points of (**) is: 


#1: x_c := x = k-^/h - s«k«-/h - t«k«-/h 

#2: y_c := y = s ■ k ■ Jh - t . k ■ Jh 

#3: z_c := z = -h + 2«s«h + 2«t«h - 4«s«t«h 

#4: y_c - x_c = (y - x = -/h<k«(2«s - 1)) 


#5: 

# 6 : 

#7: 

#S: 


SOLVE (y - x = yh-k.(2.s - 1), s) = s 
y_c + x_c = (x + y = -/h-k-CL - 2-t)) 

f 

SOLVE (x + y = ^h-k-Cl - 2-t), t) = t 


x - y - >.k ^ 
2 ■ 4 h ■ k ^ 


x + y - 7 h-k ^ 

2 ■ 4 h ■ k ^ 


EXPAND 


SUBST 


Z_c, [s, t], 


Jh ■ k 


x + y 


Jh ■ k 


2 ■ 4 h - k 


2 ■ 4 h - k 


#9: 


2 


z = 


y 


2 

k 


Derive file membran2.dfw 


2 


x 


2 

k 
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The following pictures show some buildings with hyperbolic paraboloids: 



Church St. Suitbert, Essen-Uberruhr, Architect Josef Lembroek 

http://www.tragwerkstatt.de/ipp/proiekte/kat-Q9/Q9-03.htm 




So called "Saddledome" in Calgary, Canada 


The next pictures show the restaurant LOS MANANTIALES in Xchochimilco, a suburb 
of Mexico-City. It was built by Felix Candela in 1958. The roof of this building is com- 
posed of hyperbolic parabolids. (The person on the picture is the architect of this con- 
tribution, Benno Grabinger). 


D 
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A contour plot of the hyperbolic paraboloid 


D 
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Text file for CAS-TIs: Weak Acid Speciation and Titrations 


C: conc*num_cat*cat_ch/2* (cat_ch-an_ch)->s_s_i (cone , num_cat , cat_ch , an_ch) 
C : s_s_i ( . 1 , 1 , 1 , ' 1 ) ->i 

C: 10 A ( 'a* (V(i ) / (1+V(i ) ) - . 2*i ) )^act_c_h (a , i ) 

C : pktn+a* (V(i ) / (1+V(i ) ) - . 2*i ) * (2*z-2*n+1 )->pkmn (pktn , a , i , z , n) 


C : 2 . 85^-pktl :5.7^pkt2 
C: . 51^a 

C : pkmn (pktl , a , i ,0,1 )-<-pkm1 
C : pkmn (pkt2 , a , i ,0,2)^pkm2 


C : 1 / (10 A (ph-pkm1_)+1)^mpa0(pkm1_) 
C : 1 / (10 A (pkm1_-ph)+1)^mpa1 (pkm1_) 
C : ca*mpa0(pkm1_)+ha (ca , pkm1_) 

C : ca*mpa1 (pkm1_)->am(ca , pkm1_) 




di fdata (m) 

Prgm 

Local i , d , 1_1 , 1_2 
dim(m) [ 1 ] -d 

newMat (dim(m) [1 ] ,2)^m_ 

{m[1 ,1] , m [ 2 , 1 ] , m[3 , 1 ] }->l_1 
{ m [ 1 ,2] , m [ 2 , 2 ] ,m[3,2] }-l_2 

QuadReg 1 1 , 1 2 

d ( regeq (x) , x ) | x=m[1 , 1 ] ->m_[1 ,2] 
m[1 ,1]+m_[1 ,1] 

For i , 2 , d-1 

{m[i-1 ,1] ,m[i ,1] ,m[i+1 ,1]}-1_1 
{ m [ i - 1 ,2] ,m[i ,2] ,m[i+1 ,2] }-1_2 


QuadReg 1 1 , 1 2 

d ( regeq (x) , x) | x=m [1 , 1 ]^m_[i ,2] 
m[i , 1 ]->m_[i ,1 ] 

EndFor 

{m [d-2 , 1 ] , m [d-1 ,1 ] , m [d , 1 ] }->l_1 
{m [d-2 , 2] , m [d-1 ,2] , m [d , 2] }+l_2 
QuadReg 1_1,1_2 
d ( regeq (x) , x) | x=m [d , 1 ]^m_[d , 2] 
m [ d , 1 ] -+m_ [ d , 1 ] 
m_-»fddp 

Disp "fddp calculated" 

EndPrgm 


You can find all functions together with the executable script - the first lines of the script are on top of 
the page - among the downloadable files. 
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From Subject: The Benefit of Multivalued Constants 
to Subject: Infinite Sets and Intervals 

A report about an exciting email exchange between CAS-experts 

From: Albert D. Rich 
To: David Jeffrey 

Sent: Thursday, 21 July, 2005 16:06 
Hello David, 

I whole-heartedly agree with you conclusion in your and Norman's paper entitled "Multivalued func- 
tions in computer algebra" that such functions lead to all kinds of problems and should be avoided. 

However, if properly implemented, I have been driven to the conclusion that multivalued constants are 
indispensable to good computer algebra and mathematics. For years now Derive has had multivalued 
constants for "plus-or-minus one" and for "an arbitrary point on the unit circle in the complex plane". 
For example, the equation 

ABS(z) = 1 

simplifies to z = unit_ci rcle. The attached MultivaluedConstants.dfw file illustrates that the "plus- 
or-minus one" constant is essential for doing even basic algebra. Check it out... 

Over the last several months I have incorporated the new multivalued constants #Z and #R into De- 
rive. They represent an arbitrary element of the integers and of the reals, respectively. The algebra 
required to simplify expressions involving these constants affect the entire system. Some rules are 
obvious, such as if n is an integer, 

#Z + n -> #Z 


and the not so obvious 


#Z "+ -" 1/2 -> #Z + 1/2 


and 


#Z "+ -" 1/4 -> #Z/2 + 1/4 


Of course, the real benefit of these new constant is the ability to concise and completely show all the 
solutions of equation with infinitely many solutions. For example, solving 

SIN(x) = 0 


for x used to give 


X = -7TVX = 7TVX = 0 


whereas, Derive will now give 


X = 7T*#Z 


The following are even more dramatic examples of the benefits of #Z: 


x -x 

1. S0LVE(e + e =0, x) now yields 


instead of 


X=L*(7T*#Z+7r/2) 


D-N-L#59 From Multivalued Constants to Infinite Sets, a Report p41 


3-7T- l 

V X = V X 


x x 

2. SOLVE^e = e , xj now yields 


x=1/2-V(1+8*7T*?*#Z)/2 OR x=V( 1+8*7T*?*#Z)/2+1/2 
instead of the mysterious 


i/2' w(i/(64-rr + 1) + 1) - V2) V2-WCV(64-r + 1) - 1) VMVU 

- V X = - 

4 4 


■ (7(7(64 -tt + 1) + 1) - 72) 
4 


^2 ■ L ■ -7T + 1) - 1) 

4 


j2.WWC64.r- + 1) + 1) + 42 ') 42'L'4t4t 64 -tt + 1) - 1) 

4 4 


7^- (7(7(64 -tt + 1) + 1) + ^2) 7^ 7(7(64 -tt + 1) - 1) 

+ v :< = 1 v x = 0 

4 4 


fx 2 + 1 1 

S0LVE^27 = 243 , xj now yields 


x=-7(2/3+2*7r* i *#Z/(3*LN(3))) OR x=J(2/3+2*n* i *#Z/(3*LN(3))) 
instead of 


■ G/G/(LN(3) + T7 ) - nO + -/(/(LN(3) + n ) + tt)) 1/6 ■ i.- G/G/(LN(3) + tt ) + n) - y(y(LN(3) + n ) - tt)) 

6-/CLNC3)) + 6-^(LN(3)) 

22 22 22 22 
S- (yCy(LN(3) + tt ) - tt) + y(7(LN(3) + tt ) + nO) ye- i.- (ycy(LN(3) + tt ) - tt) - 7(7(LN(3) + tt ) + tt)) 



6- 

VCLNC3)) 


- + 

6-7CLNC3)) 

ye- 

2 2 

UU(LH(3) + n ) - 

- nO + yG/(Lf 

2 2 

K3) + tt ) + tt)) 

ye- L ■ 

2 2 2 2 
^Cy(LN(3) + tt ) + tt) - y(y(LN(3) + tt ) - tt)) 


6- 

VCLNC3)) 



6.y(LH(3)) 


2 2 

cycy(LN(3) + to - 

- "■) + -/(-/(LI 

2 2 

K3) + tt ) + tt)) 

■ J 6 - L-U 

2 2 2 2 
((y(LN(3) + tt ) - tt) - y(y(LN(3) + tt ) + tt)) 


6- 

VCLN(3)) 



W(LN(3)) 


Now I am considering the ramifications of implementing multivalued constants for an arbitrary element 
of the nth roots of unity. It will be entered as #Un where n is a positive integer. This will make the 
solution of equations of the form x A n = c be the single, concise equation 

x = #Un*c A (1/n) 

instead of n messy ones. Note that as an added bonus, #U1 can replace the existing multivalued 
constant "plus-or-minus one" (of course, it will still be displayed as 

Food for thought. Look forward to your comments 

Aloha, 


Albert 
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From: Albert D. Rich 
To: David Jeffrey 

Sent: Thursday, 21 July, 2005 23:25 


Oops. Instead of #U1, below I should have said #U2 (i.e. the square-roots of unity) is equivalent to 
the "plus-or-minus one" multivalued constant. #U 1 of course equals 1 . 

Aloha 


From: David Jeffrey 
To: Albert D. Rich 

Sent: Friday, 22 July, 2005 02:09 

Hi Albert, 

At first glance, I do not see any conflict between my paper and your discoveries. The fact is that some mathema- 
tical questions do have multiple solutions. You now have the apparatus to show that to the user. This is a very 
good improvement to Derive. 

David (from hot and polluted China) 

From: David Jeffrey 
To: Albert D. Rich 

Sent: Monday, 25 July, 2005 15:52 

Hi Albert, 

Your ± example is interesting. In fact I have used that idea before. I rewrite -b ± V(b A 2-4ac) as -b ± b V( 1- 
4ac/b A 2) where the b A 2 can be pulled out exactly because of the ±. 

Some questions: if someone wants to plot the solutions, or evaluate them numerically, then you will have to split 
the symbol into two expressions. How do you think to do this ? At present Derive is happy to plot y = V(9-x A 2) 
OR y = -V(9-x A 2), so you could either transform into that format or do it behind the scenes. 

In many problems, one solves an equation and then chooses the meaningful branch. It would be good to have a 
simple way of substituting a definite sign for the ± in particular cases. 

David 

From: Albert D. Rich 
To: David Jeffrey 
Cc: David Stoutemyer 

Sent: Monday, 22 August, 2005 20:31 
Hi David, 

Thanks for your comments. 

1. Yes, I need to address the problem of plotting expressions containing multivalued constants in the 
next release of Derive. 

2. After highlighting the +- subexpression, use the Simplify > Subexpression substitution command 
(shortcut Ctrl-T) to replace it with plus or minus of the subexpression. 

Any comments about the real meat of my original email (i.e the benefits of #R, #Z, and possibly #Un)7 

Aloha, 


Albert 
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From: David Jeffrey 
To: Albert D. Rich 

Sent: Wednesday, 31 August, 2005 10:03 
Hi Albert, 

I looked again at your letter, and it seems to me there are two separate thoughts here. 

1 . An equation will sometimes have multiple solutions, and these are best expressed using your functions #R, 
#Z. Examples are evaluating an integral (solving the simplest differential equation) and solving a trig equa- 
tion. 

2. After that first step of including them in equation solutions, an equally important step is working with them. 
This is really the new part. It would be interesting to see what other systems do with #R+1, which of course 
should simplify to #R. So point 1 is I think not controversial, but point 2 is a real challenge I think. So keep 
working at it! 

David 

P.S. I may be able to drop in and visit either in December or in March. 

More later. 

From: Stoutemver, David 
To: Albert D. Rich 

Sent: Wednesday, 07 September, 2005 11:17 
Hi Albert, 

I too have come to the conclusion that we need them, but I have some uncertainties that I 
thought perhaps we could resolve by brainstorming: 

a) Do you regard, for example +-1 as the set {-1, +1} or as one arbitrary element of that set? 
If it is a set, then what about when the is written over the "+" (-+1) to mean "consistently 
take either the top or the bottom choice throughout an expression such as (a +- b)*(c-+d). 

b) How do we more generally handle correlated multi-valued constants? By this I mean 
when the same +- subexpression occurs several places in an expression. In such instances 
the "element-of interpretation yields tighter results than the "set" interpretaion. 

c) With the set interpretation, an expression such as a+-b means that we are adding a sca- 
lar "a" to a set "b", unless we regard "a" as shorthand for {a}, in which case this notational 
abbreviation might lead to apparent paradoxes. 

d) How does this kind of addition and subtraction comingle with union and set differences? 
I suspect that we shift our interpretation according to context, but I worry that being infor- 
mal and inconsistent about this will lead to difficulties. 

e) Regrettably, we never made Derive intervals first class citizens, displayable and accessible 
to the user. Do intervals denote sets or an arbitrary element of a set? How do we combine 
them with multi-valued constants? How do we combine them with Boolean expressions 
such as 0 < x <= 5? 

f) A single interval is sometimes insufficient for computing limits and representing the solu- 
tions of inequalities. The most general multi-valued constant is a "multi-interval" such as x 
= {(-3,4], 7, [9, inf)} or else x element_of {(-3,4], 7, [9, inf)}. We lose information if we gene- 
rously degrade that multi-interval into (-3, inf). 

g) With an interval, how do we indicate that we know an endpoint is tight, versus indicating 
we know that an endpoint is a bound that isn't tight, versus indicating we don't know 
whether or not an endpoint is tight? This could matter, for example, in computing limits. 

— aloha, david 
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From: Albert D. Rich 

To: David Jeffrey ; David Stoutemver 

Cc: Michel Beaudin ; Bernhard Kutzler 

Sent: Sunday, 1 1 September, 2005 23:57 

Hello David A 2, 

Thanks for your interest in the problem of multivalued expressions. We all seem to agree on the importance of 
properly using and simplifying such expressions. 

When David Stoutemyer was here this Friday, he emphasized the importance of finding a notation that 
would be natural to users. The operator is the only widely accepted way to display multivalues. 
Unfortunately, it is woefully inadequate for even some simple multivalued expressions. For example, 
an arbitrary element of the set {2, 3} in terms of is 

(5"+-"1)/2 

Representing an arbitrary element of three element sets using is even less intuitive. 

David also pointed out that most people interpret script Z and R as the set of integers and reals, re- 
spectively; not as an arbitrary element of those sets. 

To address both these problems, I propose the following: Make the lower case Greek letter epsilon 
a unary operator that reads as "an arbitrary element of". Note that the data type of the epsilon opera- 
tor is the same type as the elements of the set. For example, 

epsilon {2, 3} 

is a scalar (i.e. an arbitrary element of the set {2, 3}) that can be included in an expression like 
any other scalar. Note that epsilon of a single element set is that element. 


Since is a compact and accepted operator, internal expressions of the form 

epsilon {z, -z} 

will be displayed as "+-"z. 

Also I now propose that script Z and R denote the set of integers and reals, respectively; instead of an 
arbitrary element of those sets like I proposed before. Z and R are entered on the author line as 
#Z and #R, respectively. Then to get an arbitrary element of those sets, enter 

epsilon(#Z) and epsilon(#R) 

respectively. For example, 


SIN(z) = 0 


will simplify to 


z = 7i*epsilon(#Z) 

which reads as "z equals pi times an arbitrary element of the integers". 


Traditionally epsilon is a binary operator requiring a variable as its left operand. However, making it a 
unary operator preserves the utility of the operator, while eliminating the need to introduce a new 
dummy variable in an expression. Note that the :epsilon binary operator will continue to be used for 
domain declarations of variables. 


So in addition to finite sets (e.g. {a, b, c}), Derive now has the infininte sets #Z and #R. That still 
leaves open the question how to denote and represent arbitrary intervals (i.e. infinite sets of contigious 
points). 
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The obvious notation for denoting open, closed, open-closed, and closed-open intervals would be (a, 
b), [a, b], (a, b], and [a, b), respectively. Unfortunately, since square brackets and parentheses are 
already heavily over-loaded delimiters, using them to denote intervals could be confusing, if not am- 
bigious. For example, the notation for the closed interval [a, b] is identical to that of the two-element 
vector [a, b]. Also F times the open interval (a, b) looks like the function call F(a, b). 


Then there is the problem of how to represent intervals internally: 


An arbitrary element of the interval [-1, 1] could be represented as SIN(epsilon(#R)), and the inter- 
val itself could be represented by {SIN(epsilon(#R))}. An arbitrary element of the interval [0, 1) could 
be represented by MOD(epsilon(#R)), and the interval itself could be represented by 
{MOD(epsilon(#R))}. An arbitrary element of the interval [0, inf) could be represented by 
ABS(epsilon(#R)), and the interval itself could be represented by {ABS(epsilon(#R))}. 

Then any arbitrary closed interval [a, b] could be represented as [-1, 1]*(b-a)/2 + (a+b)/2. Any arbi- 
trary closed-open interval [a, b) could be represented as [0, 1)*(b-a) + a. Any arbitrary infinite interval 
[a, inf) could be represented as [0, inf) + a. 


Also the UNION operator can be used to represent disjoint intervals. For example, the solution of the 
inequality 

(x+ 1 )-(x — 1 )-(x — 2) > 0 


could be displayed as 


x = epsilon ((-1 , 1 ) UNION (2, inf)) 


Although cumbersome, introducing new data structures to represent arbitrary intervals directly AND 
algebraically simplifying them would be a huge project, considering the work required just to imple- 
ment the simplification of expressions involving #Z and #R. 

I hope you don't mind me using you as a sounding board. Please do not feel obligated to respond, but 
I will appreciate your comments if you do. 


Aloha, Albert 


From: David Jeffrey 
To: Albert D. Rich 

Cc: David Stoutemver ; Michel Beaudin ; Bernhard Kutzler 
Sent: Monday, 12 September, 2005 04:01 

Aloha epsilon{Albert, David}, 

I think the notation is a great idea. The standard notation x :epsilon #R becomes an abbreviation for 
x:=epsilon #R. 

One notational question that still needs addressing is that of multiple instances of the expressions. 
For example, x:= epsilon {+ 1,-1 }a - epsilon{+l,-l}b. 

Does this imply the 4 quantities a+b,a-b,-a+b,-a-b, or only the two quantities a-b, -a+b ? 

With regard to intervals, my first comment is that implementing epsilon will keep you busy for enough 
cold Hawaiian nights huddled around the fireplace. However, for reference, I checked to see what 
Maple syntax is: they have a function RealRange(a,b) which is closed unless the end point is enclosed 
in another operator Open(a). Thus the interval (a,b] is represented RealRange(Open(a),b) 

All for now: classes to teach 
David 
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From: Stoutemver, David 

To: David Jeffrey ; Albert D. Rich 

Cc: Michel Beaudin ; Bernhard Kutzler 

Sent: Monday, 12 September, 2005 07:19 

Hi Albert, 

I too like the unary epsilon idea and the generalization of union to arbitrary elements of sets. 

Regarding the user's perspective of math syntax in general and more specifically interval syntax, it is way less 
missionary work to stick as close as possible to the most commonly used notation. One way to avoid syntactic 
ambiguity with vectors and function application is to separate the endpoints of an interval with an ellipsis rather 
than a comma. For example, [3. ..5] or [3. ..5). This also emphasizes that it is a continuous range rather than 
merely the two endpoints, such as in the set notation {3, 5}. 

Regarding internal representation of an interval, Derive already has (INTERVAL ...) and (STRICT ...) analo- 
gous to Maples RealRange(...) and Open(...). This seems more direct than using sin(...), mod(...) and abs(...), 
particularly since INTERVAL is already there and STRICT gives more direct control over the endpoint type. 
Also, why not atan(...) or tanh(...) in plase of sin(...), etc.? It somehow feels wrong to have to go all the way to 
these complicated functions to represent something so fundamental as an interval. 

— keep up the good work! 

David S. 


From: Albert D. Rich 

To: David Jeffrey ; David Stoutemver 

Cc: Michel Beaudin ; Bernhard Kutzler 

Sent: Tuesday, 13 September, 2005 03:07 
Hello David A 2, 

Thanks for your comments. In response to David Jeffrey's comments and suggestions: 

1. No, the :epsilon operator is not equivalent to :=epsilon. The expression x :epsilon #R declares x 
to be real, but x will still simplify to x. On the other hand, the expression x :=epsilon #R assigns x 
the value of an arbitrary element of the reals, so x will actually simplify to the multivalued expres- 
sion epsilon #R. 

2. Multiple occurrences of multivalued constants in expressions are entirely uncorrelated (i.e. they do 
not necessarily have the same value). If you start with two independent expressions each having 
multivalued constants and then add them, why should the multivalued constants magically be- 
come correlated? 

Naturally, great care must be taken when simplifying expressions having uncorrelated multivalued 
constants. For example, if C is a multivalued expression, C*(a+b) is not equivalent to C*a + C*b (i.e. 
C cannot be distributed over or factored from a sum). More specifically, C-C is not necessarily 0! In 
general, before applying transformation rules that result in an expression being repeated, as well 
as the converse of such transformations, a check must be made to ensure that the expression is NOT 
multivalued. 

This, along with the requirements of display stepping, has forced me to move away from the philoso- 
phy of a CAS having an absolute minimal set of primitive functions. Instead a richer set is required to 
avoid the repetition of expressions. For example, TAN(z) should not be replaced internally by 
SIN(z)/COS(z) (as it is in Derive 6), since it results in two occurrences of z and is thus invalid if z is a 
multivalued expression. Similarly, 
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ASIN(z) should not be replaced by ATAN(z/SQRT(1-z A 2)), 

ABS(z) should not be replaced by z/SIGN(z), and 
MOD(z) should not be replaced by z-FLOOR(z). 

Unfortunately, each additional primitive function results in an exponential growth in the number of 
transformation rules a CAS requires. Sigh... 

In response to David Stoutemyer's comments and suggestions: 

1. Using the ellipsis along with square brackets and parentheses to designate intervals is an excel- 
lent idea! Note that Derive currently allows the ellipsis to be used to enter vectors. For example, 
[1, ..., 5] is a quick way to enter the vector [1, 2, 3, 4, 5]. Fortunately, however, the closed interval 
[1 ... 5] is distinguishable from the vector by not having commas around the ellipsis. 

2. The internal Derive utility function INTERVAL and the data structure it produces have little 
relevance to the problem of providing users with arbitrary intervals (i.e. infinite sets) that can be simpli- 
fied as first class objects. For one thing, the user can enter symbolic and/or complex endpoints for 
intervals. For example, the user can enter the interval (-#i...#i] meaning all the points on the complex 
axis between #i and -#i. INTERVAL only handles numeric and infinite endpoints. 

3. I agree that intervals are fundamental and storing them internally as {SIN(epsilon(#R))}, 
{MOD(epsilon(#R))}, etc. is too cumbersome. You propose some sort of internal interval data 
structure with slots for storing the endpoints and strictness. This would result in possibly compli- 
cated symbolic endpoints being hidden in this data structure, requiring a whole new simplification 
system. 

Instead, I propose five new constants for the intervals [0...1], (0...1), [0...1), [CL.inf), and (0...inf). Note 
that Derive already has the constant #R for the interval (-inf... inf). Then when the user enters an inter- 
val, it can be stored internally as follows: 

[a...b] stored as [0...1]*(b-a) + a 

(a...b) stored as (0...1)*(b-a) + a 

[a...b) stored as [0...1)*(b-a) + a 

(a...b] stored as [0...1)*(a-b) + b 

[a. ..inf) stored as [CL.inf) + a 

(a. ..inf) stored as (CL.inf) + a 

(-inf.. .a] stored as -[CL.inf) + a 

(-inf.. .a) stored as -(CL.inf) + a 


Then Derive's existing simplification machinery will automatically simplify expressions involving inter- 
vals. For example, the intervals 2*[a...b) and (2*b...2*a] will both simplify to the same internal form 
[0...1)*2*(b-a) + 2*a, and are thus equal intervals. 

If the leading coefficient of either of the intervals [0...1] or (0...1) is negative, the normalization rules 

-n*[0...1] ~> n*[0...1] - n 
-n*(0...1) ~> n*(0...1) - n 

are applied to make the leading coefficient positive. Therefore, [a...b] and [b...a] will both simplify to 
the same internal form, and are thus equal intervals. 
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As most recipients of this email already know, my consulting agreement with Texas Instruments In- 
corporated recently expired. Therefore whether these ideas actually make it into a future release of 
Derive is an open question. 

Anyway, thanks again for your comments. 

Aloha, 

Albert 

From: Josef Bohm 
To: Albert D. Rich 

Sent: Monday, 12 September, 2005 21:05 
Hi Albert, 

exciting exchange of mails. 

Are there any problems in publishing them in the next DNL. 

Best regards 
Josef 

From: Albert D. Rich 
To: Josef Bohm 

Sent: Monday, 13 September, 2005 10:011 
Hello Josef, 

Yes, this interaction of ideas has been an valuable learning experience for me. It brings to light how 
my thinking evolves in response to the dictates of mathematics and the comments of real world users 
like David's Stoutemyer and Jeffrey. I have no problems with you publishing them in the DNL. 

Aloha, 

Albert 


From: David Jeffrey 
To: Albert D. Rich 

Cc: David Stoutemyer ; Michel Beaudin ; Bernhard Kutzler ; Josef Bohm 
Sent: Tuesday, 13 September, 2005 12:13 

Dear All, 

1. Names. 

Strictly the symbol in question in not an epsilon. 

I have attached two expressions, printed in Knuth's Latex system, showing "x epsilon #R" and "x is-a- 
member-of #R". 

Now that Derive has Unicode, and can distinguish between "i" as a variable and "i" as sqrt(-l), it 
should also distinuish between epsilon and member-of. I am sure Knuth would disapprove of using 
epsilon. It would be the same as replacing PRODUCT and SUM with PI and SIGMA because the op- 
erators look like the letters. 
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2. Correlation. 

The expression 1/(a "+" vb) - (a - "+" vb)/(a A 2 - b) 

reduces to zero if the two "±-" operators are correlated, but not otherwise. The question is whether a 
system should offer this possibility to the users. 

In Abramowitz and Stegun, we see on page 80 the formula 

Arctan z_1 "+" Arctan z_2 = Arctan( (z_1 "+" z_2)/(z_1 z_2 "+” 1) ) 

Notice the "+" symbol. Should this be possible ? 

Michel, do you have some examples from your differential equations where the same arbitrary con- 
stant of integration appears in numerator and denominator? 

David 

From: Michel Beaudin 
To: David Jeffrey 

Cc: David Stoutemver ; Bernhard Kutzler ; Josef Bohm 
Sent: Tuesday, 13 September, 2005 23:19 

Hello my friends. 

For David Jeffrey's question to M.Beaudin: consider the separable —non linear— ODE dy/dt = (y- 
2) A 2. Then the general solution is y = 2 - l/(t+c) where c is an arbitrary constant (and y = 2 is another 
solution that does not follow from a particular choice of c). If you put the RHS of the answer over the 
same denominator, you get y = (2t+2c-l)/(t+c). Is it what you are looking for? 

Another question to All: how can we verify the formula (page 80 of Abramowith and Stegun) with De- 
rive?. In Derive, I type the following expression 

atan(x) + atan(y) - atan((x+y)/(1-x*y)). 

Simplifying gives atan(x) + atan(y) + atan((x+y)/(x*y-1 )), which is the same (atan is odd). And, plotting 
the graph (using the slide bar with y as a parameter), gives, depending on y, 2 horizontal segments. 
So, the equality seems to be valid only for a certain range. 

Best regards. 

Michel Beaudin 


From: David Jeffrey 
To: Albert D. Rich 

Cc: David Stoutemver ; Michel Beaudin ; Bernhard Kutzler ; Josef Bohm 
Sent: Tuesday, 13 September, 2005 23:35 

Michel, 

The relation is not valid for all x,y. The correct relation between the arctangents requires a step func- 
tion in the real case, or an unwinding number in the complex case. I introduced the example to show 
correlated "+ -" signs, but otherwise the example is not relevant to the discussion. 

Well in a sense it is. The lazy way to repair the relation is to add an unknown multiple of #pi, which 
could be done with Albert's new construction, but the multiple of #pi required is computable in this 
case from x,y, and that is a better solution for this case. 


David 
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From: Albert D. Rich 
To: Josef Bohm 

Sent: Thursday, 15 September, 2005 08:27 
Hello Josef, 

Thank you for your interest in publishing this exchange of emails showing how the syntax and seman- 
tics of a significant addition to Derive capabilities is arrived at through the candid exchange of ideas 
between users and developers. 

In my initial email I proposed some multivalued constants (i.e. #Z, #R, and #Un) that would greatly 
simplify the solutions returned for some equations and inequalities. Due to constructive comments 
from David A 2, this evolved into providing an arbitrary element operator for use on finite sets and the 
infinite sets script Z and R. Finally, it became apparent that primitives for the unit length intervals (i.e. 
infinite sets of reals) was also required for completeness. 

Therefore, the title of this exchange of emails should probably not be "The benefit of multivalued con- 
stants". Rather something more broad, encompassing the arbitrary element operator, infinite sets, and 
intervals. 

Aloha, 

Albert 


From: Albert D. Rich 

To: David Jeffrey ; David Stoutemver 

Cc: Michel Beaudin ; Bernhard Kutzler ; George Douros 

Sent: Thursday, 15 September, 2005 10:55 
Hello All, 

Thanks for this lively discussion on what I thought was an esoteric topic of interest to few. 

1 . I made a mistake proposing to use the lower case Greek letter epsilon to denote the arbitrary ele- 
ment (member) operator, epsilon should remain available for use as a variable name and distinct 
from the operator. Thanks to the foresight of George Douros, the Derive Unicode character set al- 
ready has a distinct character for displaying the arbitrary element operator. I plan to use ELEMENT 
as the equivalent ASCII name for the operator. 

2. Concerning the issue of correlation: DJJ gives the formula on page 80 of Abramowitz and Stegun 
as an example of the need for correlated "+-" operators. I feel the formula is just a way to show to two 
identities as one equation, thereby saving some space on the page. I have yet to see a problem for 
which a CAS would need correlation to show the solution. 

Also an arbitrary element of the set {1, -1} (i.e. "+-"1) is a simple multivalued expression for which one 
can at least contemplate correlation. However, no one would expect two occurrences of arbitrary 
elements of more complicated sets (e.g. the set of all reals or intervals) to be correlated. If correlation 
is required, use a variable with an appropriately declared domain. For example, the familiar identity 

SIN(x) A 2 + COS(x) A 2 = 1 

would be confusing if written as 

SIN(ELEMENT(#R)) A 2 + COS(ELEMENT(#R)) A 2 = 1 
with the assumption that the ELEMENT(#R)s are correlated. 
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3. The 5 primitive intervals proposed in my last email can and should be reduced to just the 3 unit 
length intervals by eliminating the 2 infinite length intervals [0...inf) and (0...inf). Then intervals of the 
form [a. ..inf) can be replaced with [0...1 )*inf + a, and intervals of the form (a. ..inf) can be replaced with 
(0...1 )*inf + a. 

Aloha, 

Albert 

From: David Jeffrey 
To: Albert D. Rich 

Cc: Michel Beaudin ; Bernhard Kutzler ; David Stoutemyer ; George Douros 
Sent: Thursday, 15 September, 2005 14:49 


With regard to 2. in Albert's mail: 

In Michel's example of an ODE whose solution had the same constant in numerator and denominator, it ap- 
peared that the constant only appeared in 2 places because Michel put the answer over a common denominator. 

Hence it looks as though a CAS should be able to get every answer in a form in which the constant only appears 
once. However, I am sure that if Michel looks around a little longer, he can find us an example in which it is not 
possible to put the answer in such a form. Even if he cannot find an example right away, there remains the pos- 
sibility that a CAS would also put the expression over a common denominator. Perhaps the best way to handle 
correlation if it is needed is to introduce an auxiliary variable. To create a quick example: 


simplifies at present to 


S0LUTI0NS([x + y = 1], [x, y]) 

[[@ 2 , 1 - @ 2 ] ] 


where at present Derive does not state @2 : ELEMENT #R. If it did, however, then all occurrences of 
@2 would have the same. If however this were replaced with 

[ [ELEMENT(#R) , 1- ELEMENT(#R) ] ] 


then the two objects would be uncorrelated. 
DJJ 


From: Albert D. Rich 

To: David Jeffrey ; David Stoutemyer 

Cc: Michel Beaudin ; Bernhard Kutzler ; George Douros 

Sent: Sunday, 25 September, 2005 07:57 

Subject: Infinite sets and intervals 
Hello All, 

At the end of my email dated 12 September 2005 below, I propose representing general inter- 
vals internally using 3 primitive unit interval constants multiplied by a scale factor and added to by an 
offset. The alternative is to have an interval function having as arguments the bounds and strictness. 
The reason I give for using unit interval constants was that "Derive's existing simplification machinery 
will automatically simplify expressions involving intervals." 


p52 


A Call for Ideas 


D-N-L#59 


Now I have more rationale for my belief: The complex number x+y*#i can be represented by a imagi- 
nary unit constant #i multiplied by y and added to by x. The alternative is to have a REIM function 
having as arguments x and y. Unfortunately, using REIM forces a CAS to properly handle two forms 
of sums and two forms of products all over the place! For this reason, David Stoutemyer and I have 
been driven to the conclusion that using REIM instead of an imaginary unit constant is significantly 
more difficult to implement. 


The choice between using unit interval constants instead of an interval function is entirely analogous 
(isomorphic?) to using an imaginary unit constant instead of a REIM function. Thus unit interval con- 
stants should be used, qed 


Aloha, Albert 


P.S. I took the liberty of changing the subject of this email exchange since the idea of "multivalued 
constants" has evolved into a much more ambitious plan to introduce an arbitrary element operator 
and support for infinite sets and intervals. 


Call for Ideas for a DERIVE - Competition 

We'd appreciate ideas and comments from DUG members to help us devise a Derive competition. 

The competition has two aims: 

1. To raise awareness of the benefits of using Derive amongst newcomers. 

2. To encourage existing Derive users to make better educational use of Derive. 

Therefore the competition will have two questions. 

The question for newcomers should be one which could be answered by someone using the 30-day 
trial version, so we might stimulate new users to try Derive. Perhaps it might be something like "Cre- 
ate a Derive worksheet showing a graphical and algebraic approach to solving simultaneous equa- 
tions." 

The prize for this might be a free Derive single user licence. 

The question for existing users, which could win them a School Site Licence Plus (or equivalent value 
licence) could be something like one of the following. 

Project: Meaningful Use of Slider Bars (2D & 3D) 

Project: Calculate in the Algebra Window and present in the 3D- or 2D-Plot-Window 

Project: Combine CAS & Programming 

Project: DERIVE & Voyage 200 - A Game for Two 

Or perhaps a mathematical challenge. 

Do any DUG members have suggestions for questions which would best achieve these aims? For the 
widest impact perhaps the level should be school level. 

Cheers 
Philip Yorke 

Chartwell-Yorke Ltd. 

1 14 High Street, Belmont Village, Bolton, Lancashire, BL7 8AL 
England, tel (+44) (0)1204 811001, fax (+44) (0)1204 811008 
info@chartwellyorke.com, www.chartwellyorke.com 


